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Abstract 



We develop an analytical theory to explain the experimentally-observed morphological transitions 
of giant vesicles induced by AC electric fields ([1]). The model treats the inner and suspending media 
as lossy dielectrics, and the membrane as an ion-impermeable flexible incompressible- fluid sheet. 
The vesicle shape is obtained by balancing electric, hydrodynamic, and bending stresses exerted 
on the membrane. Considering a nearly spherical vesicle, the solution to the electrohydrodynamic 
problem is obtained as a regular perturbation expansion in the excess area. 

The theory predicts that stationary vesicle deformation depends on field frequency and conduc- 
tivity conditions. If the inner fluid is more conducting than the suspending medium, the vesicle 
always adopts a prolate shape. In the opposite case, the vesicle undergoes a transition from a 
prolate to oblate ellipsoid at a critical frequency, which the theory identifies with the inverse mem- 
brane charging time. At frequencies higher than the inverse Maxwell- Wagner polarization time, 
the electrohydrodynamic stresses become too small to alter the vesicle's quasi-spherical rest shape. 
The analysis shows that the evolution towards the stationary vesicle shape strongly depends on 
membrane properties such as viscosity. The model can be used to rationalize the transient and 
steady deformation of biological cells in electric fields. 
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1 Introduction 

Electric fields are widely used for cell manipulation. Weak fields influence cell signaling, wound 
healing, and cell growth ()2l |3l HI |5]) . Strong pulsed fields can induce transient perforation of the 
cell membrane, which enables the delivery of exogenous molecules (drugs, proteins, and plasmids) 
into living cells (j5l [6| . 

Biological cells exhibit various frequency-dependent behaviors in AC electric fields: orientation, 
translation (dielectrophoresis), and rotation. These phenomena have stimulated considerable mod- 
eling effort aimed at understanding of the physical mechanisms of the interaction of electric fields 
with cells and tissues. A common theme among different theoretical models is the assumption that 
the cell is a lossy dielectric particle of /ixed shape (a sphere (f71 l8l l9l [TO\i or an ellipsoid fjll^ il2,) ). For 
example, the orientation of cells can be predicted by considering the torque on an ellipsoid due to 
the effective dipole moment induced by the electric field (|T3l [TS|) ; the dipole based theory has been 
successfully applied to interpret electro-orientation of erythrocytes (fT5|l . 

Cells, however, are soft objects, which deform when subjected to electric fields. The cell mem- 
brane plays a critical role in this process. A number of studies have focused on the membrane 
shear elasticity because of the interest in the mechanics of the red blood cell (|16| [T71 [T8|) . The 
lipid bilayer is the main structural component of the cell membrane, yet, the electrodeformation 
of closed pure lipid bilayer membranes (vesicles) has been considered only to a limited extent 
(jl9t [501 EH)- There is increasing interest in this problem, particularly in relation to electroperme- 
abilization (|22| l23l [2^ l25l l26l 1271 ESI- Recent experiments have shown that vesicle behavior in 
electric fields exhibits peculiar features. Vesicles subjected to a direct-current-electric pulse can 
deform into elliptical (I29p or cylindrical shapes ()30p . Alternating-current-electric fields deform 
vesicles into prolate or oblate ellipsoids depending on the frequency and the conductivities of the 
interior and suspending fluids ([H EH l32ll . 

The physical mechanisms responsible for vesicle electro deformations are not fully understood 
at present time. A prolate shape can be explained by the electric pressure pulling the vesicle at 
the poles. However, the oblate shapes remain an open problem, in particular, the fact that they 
are observed only when the conductivity ratio of inner and outer fluids is less than unity (152]) . 
Attempts to explain the oblate shapes have been made (I33p . Peterlin et al. (I34p showed that 
anisotropy in the dielectric constant of the membrane could lead to oblate shapes, but this model 
does not account for the observed dependence on the conductivity ratio. Hyuga et al. (J35J) realized 
that the fluid environment is not just a passive milieu and that the dynamic coupling between 
changes in membrane shape and motion in the surrounding fluids is important. However, the fluid 
flow in their model was described by an ad hoc equation that includes fluid acceleration. On the 
contrary, fluid motion at the micro-scale is in the viscous Stokes regime, where friction effects 
dampen fluid acceleration ()36p . Thus, the validity of their analysis is questionable. 

In this paper we develop a model for vesicle dynamics in electric fields. The transient vesicle 
deformation is determined by evaluating the forces exerted on the membrane (|37p . Our theory 
builds on the large amount of research devoted to the electrohydrodynamics of drops ([551 1321 [^0]) . 
Drops can adopt oblate shapes because the electric field causes continuous fluid flow, which pushes 
fluid toward the equator. The physical mechanism behind the electrohydrodynamic flow is the 
following: If the fluids are leaky dielectrics, i.e., possess finite conductivity, free charges accumulate 
at the drop interface. The electric field acting on these changes creates a tangential electric force, 
which drags the fluids into motion (see Figure 1). 

Notwithstanding the qualitative similarity between drops and vesicles, the extension of the 
"leaky dielectric" model from drops to vesicles is not a straightforward task because the mechanics 
of lipid membranes is far more complex than the mechanics of fluid- fluid interfaces. There are two 
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major challenges. First, the lipid membrane is essentially an insulating shell impermeable to ions. 
When an electric field is applied, charges accumulate on both sides of the bilayer and the vesicle 
acts as a charging capacitor. Second, since the lipid bilayer contains a fixed number of molecules 
the membrane is incompressible. Under stress, the membrane develops tension, which adapts itself 
to the forces exerted on the membrane in order to keep the local area constant. At steady state 
the gradients in tension counteract the tangential electric force and the electrohydrodynamic fiow 
stops. Our approach rigorously accounts for these phenomena. 

The paper is organized as follows: Section [2] describes the physical model and formulates the 
governing equations. Section |3] outlines the solution and discusses the frequency dependence of the 
electric stresses, and Section |4] shows the theory predictions for the vesicle shape as a function 
of frequency, conductivity ratio and other physical parameters of the system such as membrane 
viscosity. 



2 The model 

2.1 The physical picture and characteristic time scales 

Let us consider a vesicle with no net charge formed by a membrane with conductivity Amm) dielectric 
constant emnu and surface viscosity r]mm- The bilayer thickness is about h ~ 5nm, thus on the 
length scale of a cell-size vesicle (radius a ~ lO/um) the bilayer membrane can be regarded as a 
two-dimensional surface with capacitance Cm = ^mm/h and conductivity Gm = ^mm/h. The vesicle 
is filled with a fiuid of viscosity rji^, conductivity Ajn, and dielectric constant ein, and suspended in 
a different fiuid characterized by ryex, Aex, and eex- The physical properties of the fiuids and the 
membrane are assumed to be frequency-independent. 

The vesicle is subjected to a uniform AC electric field with an amplitude £"0 

= Eocos{iot), (1) 

where oj is the angular field frequency and t is the time. Free charges move under the action of the 
electric field. The density of free charges in the bulk fiuids decays on a time scale (|39| HOl) 



tc = ^ , i = in,ex (2) 

where e and A denote the absolute permittivity and conductivity of the fiuid. Hence, for frequencies 
Lu < free charges are present only at boundaries that separate media with different electric 
properties. The rate of accumulation of charges at the interface of a macroscopic object, e.g., a 
sphere, is given by the Maxwell- Wagner polarization time (|lUp 

, ^in + 2eex /„v 

tMW = T , n,^ ■ l-^j 

Ain -t- /Aex 

The electric field acts on the free charges at the interface and gives rise to a force, which is 
tangential to the interface. In the case of a simple fiuid-fiuid interface, e.g., a drop, only a hydro- 
dynamic force can balance the shearing electric force. As a result, the fiuids are set in continuous 
motion, the so called electrohydrodynamic (EHD) flow (|38p . The EHD flow is characterized by a 
time scale, which corresponds to the inverse of the shear rate imposed by the tangential electric 
stress 

^el = • (4) 
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In the case of drops, the flow inside is toroidal with a direction either from or towards the poles 
depending on the surface charge distribution as illustrated in Figure [T} 

The membrane represents a more complex boundary compared to fluid- fluid interfaces. First, 
it is an insulating shell and charges accumulate on both the inner and outer surfaces, as illustrated 
in Figure [2} Hence, a vesicle of radius a acts as spherical capacitor that charges on a time scale 
given by dU USD 

tcap = ClCm I T h TTT — ) • (5) 

Second, the membrane mechanics is governed by resistance to bending. A distortion of the mem- 
brane shape relaxes on a time scale 

= , (6) 

Hi, 

where k is the bending modulus. The curvature relaxation depends on the viscosity of the surround- 
ing fluids, because the membrane has to move fluid in order to return to its preferred configuration. 

It is instructive to estimate the magnitude of the characteristic time scales involved in vesicle 
electro deformation. Typical experimental conditions involve fluids with conductivities in the range 
A ~ 10-3 - and electric fields of the order of E ~ IkV/cm (jU |M1 EOl E21 IMD- In 

physiological conditions, e.g., blood, the internal conductivity of an erythrocyte is much higher 
~ 0.5S/m (|15p . The typical size of a vesicle or cell is a ~ 10//m. The inner and outer fluids are 
essentially water: viscosity r] ~ 10~^Pa.s,and density p ~ lOOOkg/m^. The membrane capacitance 
is Cm ~ W~'^F/m? (143p and bending rigidity k ~ W~^^J. Therefore, for vesicles, we estimate the 
basic charging time and the Maxwell- Wagner polarization time ~ tMW ~ the membrane 

charging time tcap ^ lO^^s, the electrohydrodynamic time td ^ lO^^s, and the bending relaxation 
time tK ~ 10s. 

We see that vesicle dynamics in electric fields involves processes that occur on very different 
times scales. Vesicle deformation takes place concurrently with fluid motion. The electric field 
adjusts to a new boundary configuration much faster than the fluid moves, because conduction 
(and hence charge redistribution) is fast, Imw ^ ^ei- Hence, the electric field depends only on the 
instantaneous vesicle shape; it is quasi-static. The flow time scale is comparable to the capacitor 
charging time, igi ~ tcap- The interplay between these two time scales is responsible for the observed 
dynamics of vesicles in electric fields ([1]) . 



2.2 Governing equations 

In essence, the leaky dielectric model consists of conservation of current, which obeys Ohm's law, 
and the Stokes equations to describe fluid motion (j4tJ) . Charges carried by conduction accumulate 
at interfaces, and bulk phases become charge- free on a very fast time scale given by tc, Eq. |2} 
Accordingly, electromechanical coupling occurs only at boundaries. 



2.2.1 Electrohydrodynamic problem 

Electric field: In the absence of bulk charges, the electric potential, ^, for a quasi-static elec- 
tromagnetic field is the solution of the Laplace equation 

V2^> = 0, E = -V^>. (7) 
The membrane acts as a capacitor. Accordingly, the potential undergoes a jump across the interface 

$in_$cx^^^(^^^^ at r = rs, (8) 
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where r = rg denotes the position of the interface in a coordinate system centered in the vesicle 
(see Fig. |2]). The relation between the transmembrane potential and the membrane capacitance 
depends on geometry. The spherical shell is a widely used model for cells and vesicles (j9l E]) . 
although a spheroidal geometry has also been considered (llip . Free charges at the interface cause 
discontinuity in the normal component of the displacement vector 

n • [eexE^" - einE''^] = Q{uj, t) at r = r„ (9) 

where n is the outward unit normal vector and Q is the free charge density. Neglecting effects of 
charge convection along the surface by fluid motion, the conservation of electric currents at the 
interface requires that 

n • [AexE^^ - AinE'-] = at r = r,. (10) 

The forces due to an electric field E are calculated from the Maxwell stress tensor 

T^^ = e (EE - -E^I ) , (11) 
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where I denotes the unit tensor. A harmonic electric field can be written as 

^cos(cut) = ^ (E + E*) . (12) 

It gives rise to a non-zero time-averaged component of the Maxwell stress tensor 

Tf{uj) = ^ [EE* + E*E - I^I^I] , (13) 

which is responsible for the steady deformation of the vesicle. 

All electric variables (electric field, potential, charge density) vary harmonically with time 
u{r,t,u;) = u(r, a;) exp(ia't). Hence, hereafter unless specifically stated, we will always refer to 
the amplitude of an electric variable, u(r,u;), and we will omit the bar for convenience. 

Hydrodynamic field: Vesicle deformation is accompanied by motion in the surrounding fluids. 
The fluid velocity, v, and pressure, p, inside and outside the vesicle are described by the Stokes 
equations ([Ml Si} 

^ ^ y . rphd V • V = , (14) 

dt ^ ' 

where the bulk hydrodynamic stress is 

T^'i = -pi + r?[Vv + (Vv)t] , (15) 

where the superscript f denotes transpose. 



Eq. 14 is a simplified version of the more general Navier-Stokes equation. First, inertial effects 
are neglected because at the length-scale of the cell water is effectively very viscous. Second, 
the bulk stress has no contribution from the electric field because there are no free bulk charges. 
Moreover, the unsteady term d^r /dt can be neglected if the diffusion of momentum, t„ = a?p/rj^ is 
faster that the changes in the electric field, i.e., u) <t~^ (HS]) . The linearity and quasi-steadiness of 
the Stokes equations, and the decoupling of the electric and hydrodynamic equations in the bulk 
greatly simplify the solution of the problem. 
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Far away from the vesicle, the fluid is at rest and the How field vanishes, v*^^ 0. Velocity is 
continuous across the interface 

^in^^ex^^^ at r = rs. (16) 
The interface moves with the fluid (I46p 

§^=v,.„. (17) 

Electromechanical coupling: The vesicle shape is determined by the balance between electric, 
hydrodynamic, and membrane tractions (surface force densities) at the interface r = Vg 

n • [(??exT^^^'^^ - m^T'^'^n + (ecxT^''^^ - einT^i'''^)] = r'"^'" , (18) 

where flexoelectric bending of the lipid bilayer is neglected (|47| HSjl . For example, at rest, when the 



electric field is off, Eq. 18 reduces to the familiar Euler-Lagrange equation (149 p . which states that 



there can be a jump in the hydrostatic pressure across a membrane due to membrane tractions 

pin _ pCx ^ 2aH - K [AH^ - AKH + 2V^iJ] . (19) 

where k is the bending rigidity, H and K are the mean and Gaussian curvatures. In the next 
section we discuss the membrane stresses in more detail. 



2.2.2 Membrane mechanics 

The pure lipid membrane consists of two sheets of lipid molecules. The molecular thickness imparts 
resistance to bending. Within the framework of the minimal model (.50), the bending resistance 
gives rise to a surface force density 

r'^ = -K [AH^ - AKH + 2V^i7) n . (20) 

The surface gradient operator is defined as Vg = Is • V, where the matrix Ig = I — nn represents a 
surface projection. 

The membrane leaflets consist of fixed number of lipids, which are optimally packed with fixed 
area per lipid (under moderate stresses). As a result, a membrane element only deforms but can 
not change its area. Under stress, the membrane develops tension, which adapts itself to the forces 
exerted on the membrane in order to keep the local and total area constant. Hence, the tension 
is non-uniform along the interface and varies with forcing. The membrane tension gives rise to 
surface force density 

t" = 2aHTi-Vga . (21) 

where a denotes the local membrane tension 

Lipid molecules are free to move within the monolayer, and therefore, in contrast to solid-like 
polymerized membranes, the lipid bilayer membrane is fluid with a zero shear-elastic modulus (T3T]). 
How fast the membrane flows, however, depends on the rate-of-deformation. The viscous stresses 
developing in the flowing membrane are 

r" = v^rA^s^s + (V,v,)t] . (22) 

The membrane viscosity of lipid bilayers is relatively low, r/mm ~ 10~^A^s/m, and its effects are 
usually negligible. However, polymersomes can have very viscous membranes characterized by large 
membrane viscosity Tymm ~ 10~^A^s/m {iST \ IH ^ . 
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2.3 Dimensionless parameters 

It is more convenient to describe the problem in terms of non-dimensional parameters. Casting 
equations in dimensionless form helps show the generality of application to a broad class of situations 
rather than just one set of dimensional parameters. 

Henceforth, bending stresses and tension are normalized by /t/a^; all other quantities are 
rescaled using t^cx, Cexi Aex, a, and Eq. The fluid velocity scale is vq = eex-E'o^/??ex- The elec- 
tric and viscous stresses are rescaled by Cex-E'o- Time and frequency are nondimensionaized with 
the basic charging time tc = Cex/Aex- 

The electric capillary number compares the shape-preserving bending stresses to the shape- 
distorting electric stresses, 

Ca=^ = ^^?^^. (23) 
The other relevant parameters are the ratios of the electric properties of inner and outer fluid 

A = ^, S = ^ (24) 

'^ex £ex 

and the viscosity ratio 

X = — • (25 

Vex 

The dimensionless membrane conductivity and capacitance per unit area are 

^ _ Amm ^ _ emm 

<-^m — — T — ) <-^m — ) i^Oj 

where the dimensionless membrane thickness is x = h/a. The importance of membrane viscous 
stresses is reflected by the magnitude of the surface viscosity parameter 

X. = ^ . (27) 

The surface viscosity parameter for a lipid vesicle is relatively small, Xs ~ 1, but for polymer- 
somes can be quite large, Xs ~ 10'^. 

We estimate that Ca ~ 10^ 3> 1 from the typical values discussed at the end of Section |2]l. The 
dielectric constants and viscosity ratios are 5, x ~ 1 • In physiological environments and biological 
applications, the conductivity ratio can vary between 10~^ and 100. 



3 Solution for small deformations 

In a coordinate system centered at the vesicle, the radial position of the vesicle interface is 
described by 

r, = l + /(0,0), (28) 

where / is the deviation of vesicle shape from a sphere. For a nearly spherical vesicle, f <^ 1. In 
this case, the exact position of the interface is replaced by the surface of a sphere of equivalent 
volume, and all quantities that are to be evaluated at the interface of the deformed vesicle are 
approximated using a Taylor series expansion. The solution for electric and flow fields is derived 
as a regular perturbation expansion in some small parameter, e.g., the excess area. 

In this study we perform the leading order analysis. Accordingly, the electric and hydrodynamic 
fields are evaluated about a sphere. First, we determine the electric field and the electric tractions 
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(surface force density) exerted on the membrane. Second, we determine the hydrodynamic tractions 
needed to satisfy the force balance Eq. [18] and the corresponding velocity field. Finally, we use the 
kinematic condition Eq. [17] to find the shape evolution. 



In Eq. 28 the function / representing the perturbation of the vesicle shape depends only on 
angular coordinates. Thus, it is expanded into series of scalar spherical harmonics Yjm given by 



Eq. 47 in Appendix [B] 

oo j 

f ~ ^ ^ ^ ^ fjmYjm 1 (29) 
j=2 m=-j 

Solutions for the electric field are growing and decaying harmonics which derive from V(r'^Tj^) 
and V{r~^~^Yjm)- The uniform applied electric field along the z-direction, defined by Eq. hi is 
described by the j = 1 harmonic 




V 3 

Accordingly, the induced electric field has j = 1 , m = symmetry. 
3.1 Electrostatic field and stresses 

The model for the electric field is based on the classic works by Schwan and coworkers ()4ip . They 
have shown that an external AC electric field induces a potential across the membrane of a spherical 
sheh dSSI, 

A^> = y„(u;)cos6', (31) 

where 

3 1 

^"^^""^ = 2 1 + (G„ + iu;C™)(l/A + 2) " ^^^^ 

The transmembrane potential is very sensitive to the membrane thickness. Figure [3]a illustrates 
the variation of the transmembrane potential with frequency for a vesicle with a fixed size and 
two values of the membrane thickness, corresponding to a giant unilamelar lipid vesicle and a 
polymersome. For a simple fluid-fluid interface (a drop, i.e., x = 0), the transmembrane potential 
is zero. 

The electric tractions exerted on the membrane have radial and tangential components 

T^' = rf [1 + 3 cos(2e)] f + Tf sm{2e)e . (33) 
In terms of the electric field, the electric pressure can be written as 

rf = \[{El-f - {ETf - S{{E^f - {E^f)] , (34) 

and the tangential electric force is 

rf = EtQ + SE^Vm{u:) sin 9 , (35) 

where we have used the definition of surface charge Q Eq. ]9j The amplitudes of the electric 
tractions, rf- and r^', depend only on the electric properties of the media. Their expressions are 



given by Eq. p2]and Eq. 53 in Appendix C.l 



The electric stresses are complicated functions of the frequency oj as illustrated in Figures ]4ja 
and]4]b. We can distinguish three regimes: 



9 



Low frequencies, uj < wi; In this case, the membrane shields the vesicle interior and the 
electric field inside zero, as seen from Figure [Sjb. The electric pressure is positive at the poles, 
and negative at the equator, thus pulling the vesicle into a prolate shape. The tangential electric 
stress is zero everywhere on the surface, because both induced charge and internal electric field are 
zero. In contrast, the tangential electric stress at a simple fluid-fluid interface, i.e., zero-membrane- 



thickness, is non-zero even at low frequencies, see Eq. 55 The electric pressure changes sign and 



the tangential electric traction becomes significant above a frequency wi given by (j9l HT]l 

Gffi 2i\. (3(3) 

Cm Cm{-^ + 2) 

which reduces to 1/tcap, Eq- [S] if the membrane is non-conducting. 

Intermediate frequencies, lvi < < UJ2: In this frequency window, the membrane capacitor 
becomes "short-circuited" and the vesicle interior participates in the conduction process. The 
onset of decrease in the transmembrane potential and increase in the interior electric field coincides 
with the appearance of tangential electric tractions and negative electric pressure, as seen in Figures 
[3]and|4j The tangential electric stress is mainly due to the free charges on the membrane (|38p , see 
Eq. [35J Because of the different conductivities of the inner and outer fluids, charges accumulate 
at different rates on the membrane physical surfaces. Charge densities on the inner and outer 
membrane surfaces can become imbalanced, which gives rise to a non-zero effective interfacial 
charge density as shown in Figure [5] The effective charge is zero at low frequencies because the 
membrane capacitor is fully charged, having equal charge densities on the inner and outer membrane 
surface, and at high frequencies because of insufficient time for interface charging. 

High frequencies, uj > lu2- The inverse Maxwell- Wagner polarization time, Eq. [3j defines a 
critical frequency 

(3.) 

above which tangential stress starts to decrease. It vanishes at very high frequencies, where all 
media behave as perfect dielectrics. In this frequency regime, the electric pressure is small, but 
positive with magnitude ~ (5 — 1)^/(5' -|- 2)^, which leads to small prolate deformation. 

3.2 Hydrodynamic field and vesicle deformation 



The stress balance at the interface Eq. 18 shows that the electric tractions need to be compensated 
by membrane and hydrodynamic forces. The latter can be found using the general solution for a 
nearly spherical vesicle subject to an external field of arbitrary symmetry developed in Ref. (I5l]) . 
Details of the solution are presented in the Supplementary material. 

The vesicle area. A, exceeds the area needed to enclose the volume of the interior fluid, Aira^. 
At rest, the excess area is redistributed among all shape modes 

A = A/a^ -47r = ^Yl -^^^ " ^^^^ + • (38) 

j=2 m=-j 

Therefore, in order to accurately describe vesicle deformation, in general, we need the evolution 
equations for all shape modes. These are derived in (|54p (see also the Supplementary material) and 
have the general form 

^ = C,m + Ca-\Ti + aoT2)fjm + O (f^) . (39) 
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The first term describes the distortion of the vesicle shape by the electrohydrodynamic flow. The 
term including Ca is associated with shape relaxation driven by the membrane stresses. The 
coefficients Cjm, Fi, and T2 are listed in the Supplementary material. The effective tension (Tq 
depends on the vesicle shape, which in turns depends on the applied electric field. 

In general, the apparent area of a vesicle, is lower than its true area. A, because of suboptical 
fluctuations in the shape modes. For example, a quasi-spherical vesicle at equilibrium is character- 
ized by a zero apparent area, i.e., A = 0. However, even though the membrane in inextensible, the 
vesicle can deform and increase its apparent area due to flattening of the shape undulations. This 
leads to an increase in the the membrane tension (15511 



o-Q = (Jin exp 



1 



(40) 



where Um is the initial membrane tension. 

In the next section we simplify the theory for the case of vesicle electrodeformation induced by 
an uniform AC electric field. 



4 Results 



4.1 Deformation of a quasi-spherical vesicle 

When the electric field is turned on, it generates electrohydrodynamic flow with the same symme- 
try as the electric stresses. The corresponding fluid velocity, which is responsible for the vesicle 
deformation, is given by 

' C20 = 8,/; .. (41) 



5 23x + IQXs + 32 



where the electric stresses are given by Eq. [52] and Eq. 53 in Appendix C.l, and the viscosity 
parameters x and Xs are defined by Eq. |25] and Eq. [27] Since electric stresses directly affect only 
the ellipsoidal j = 2, m = mode, the most important contribution to the vesicle deformation 
comes from the "elongational" /20 mode. Moreover, because the shape modes are coupled through 
the area constraint Eq. |38j the area stored in the j 7^ 2 modes is transferred into the ellipsoidal /20 
mode. The maximum possible vesicle deformation corresponds to elongation where all excess area 
is stored in the /20 mode 

/2T" = ±\/f> (42) 



where a positive sign corresponds to a prolate deformation. 

The shape evolution strongly depends on the effective tension ctq. For a quasi-spherical vesicle. 



using the relation between excess area and shape modes Eq. 38 and including only the dominant 
contribution from the /20 mode, we can rewrite Eq. [40] as 

4k 



o-Q = CTm exp 



/20 



(43) 



Inserting into Eq. 39 we obtain that the shape evolution of a vesicle in AC electric field is described 
by the following non-linear equation 



20 



- Ca 



-1 



24 



6 + cJi„exp (^T^floit) 



/2o(t) 



(44) 



dt " 23x + 16xs + 32 

Our theory can also be applied to vesicles with non-spherical rest shapes, i.e., non-zero initial 
apparent shown in Appendix IdI 
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4.2 Discussion 



The shape evolution obtained from Eq. 44 for several frequencies is illustrated in Figure [6}a. 



The vesicle deforms on a hydrodynamic time scale approximately given by = The time 

needed to reach stationary shape depends strongly on the viscosity contrast between the inner 
and outer fluid. Figure [6|b shows that increasing the viscosity of the inner fluid slows down 
the shape evolution. The viscosity effect may become important in the electrodeformation of 
red blood cells, which are characterized by x 10. Another factor that can slow down the 
shape evolution even more dramatically is the membrane viscosity, as illustrated in Figure |6jc. 
In the case of polymersomes or lipid membranes undergoing fluid-to-gel transition, the membrane 
viscosity parameter can reach values of the order of 100 (I56p . The sensitivity of shape evolution to 
membrane viscosity suggests a novel method for determination of the membrane viscosity where 
the experimental effort is minimal. 

The steady shape of a vesicle in AC electric field is calculated by evaluating Eq. [44} Figure [7] 
illustrates the steady shapes of vesicles in AC field as a function of frequency for different conduc- 
tivity ratios. The theory predicts that the type of deformation, prolate or oblate, is determined 
primarily by the frequency and the conductivity ratio. At low frequencies to < ivi the deformation 
is prolate. For frequencies lv > oJi vesicles are prolate or oblate depending on the conductivity 
ratio. At even higher frequencies, the deformation becomes again prolate but very small and the 
vesicle appears spherical. This observation corresponds well to the experimental data. 

Next we analyze these morphological transitions in more detail. 

4.2.1 Prolate-oblate transition for A < 1 at low frequencies 



The transition frequency cui corresponds to the capacitor charging time Eq. 36 



At low frequencies, u; < wi, vesicle deformation is due solely to the positive electric pressure. It 



is maximal at the poles, see Eq. 52 in Appendix C.l The vesicle is pulled apart and thus adopts 
a prolate ellipsoidal shape. 

At uj > u>i, the tangential electric traction becomes significant and the electric pressure is 
negative, as seen from Figure |4j The shearing tangential force induces electrohydrodynamic flow, 
similar to the one observed with drops (Figure [T]) . If A/S" < 1 the flow is directed from the poles to 
the equator and the resulting deformation is oblate; if A/S > 1 the flow is directed from the equator 
to the poles and the resulting deformation is prolate. Therefore, prolate-oblate transition is possible 
only if A/5 < 1. In experiments with vesicles ([II [30 l [32 |) . the inner and outer fluids are sucrose and 
glucose, which have similar dielectric constant, S" ~ 1. Oblate shapes were reported for conductivity 
ratio less than 1, in agreement with the condition A/S < 1. In the case of biological cells, the 
difference between the dielectric constants of the cytosol and the cell environment is also small, and 
therefore similar deformation behavior is expected. In the case of drops, the electrohydrodynamic 
flow persists for as long as the electric field is applied because only viscous stresses can balance the 
tangential electric surface force. In contrast to drops, the electrohydrodynamic flow in vesicles is 
not sustained. It stops when the vesicle reaches steady deformation because the membrane tension 
counteracts the electric tangential force. 

The capacitor charging time decreases with the size of the vesicle. Therefore, the smaller the 
vesicle, the higher the transition frequency. For nanometer size vesicle this frequency is in the MHz 
range. Thus, nano-vesicles are expected to deform only into prolate ellipsoids when subjected to 
AC fields with frequency less than a MHz or DC pulses with length longer than 1 ^s, which is in 
agreement with experimental observations (|57p . The theoretical predictions for the prolate-oblate 
transition frequency observed for giant vesicles are in good agreement with experiments ([T]), as 
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shown in Figure [8j Note that the reported experimental data was collected for vesicles with various 
sizes and conductivity conditions. 

4.2.2 Oblate-prolate transition for A < 1 at high frequencies 

At high frequencies to ^ toi, the transmembrane potential vanishes, as seen in Figure [3} The 
electric tractions are given by the zero-thickness results Eq. [54] and Eq. 55 in Appendix |C.1[ The 



forcing term in the shape evolution Eq. 44 changes sign at a frequency 

"45-(A + l)2n ' 



no 



(S-l) 



(45) 



Correspondingly, the vesicle deformation changes from oblate to prolate at this frequency. The 
transition frequency becomes very large when the dielectric constants of the fluids are compa- 
rable. For vesicles filled with sucrose and suspended in glucose solutions this frequency is about 
lOMHz, which is in the frequency range where electric tractions have already become too small 
to deform the vesicle. Thus this oblate-prolate transitions was not observed in the experiments of 
Aranda et al. (jT|; instead, the vesicles remain spherical. Thus far, the prolate-oblate transition 
has been reported only for drops 

If 0,2 < ^1, the oblate deformation would be impossible. This situation arises if the membrane 
becomes highly conducting, e.g., because of poration. Another possibility is a thick membrane or 
small vesicle with 

h {A + 2)^Sm , . 

a {S + 2){2A + Gm{A + 2)) ^ ' 

where Sm = emm/cex- For a typical bilayer thickness of 5nm this condition holds for vesicle size 
below lOOnm. This prediction is in agreement with experimental studies of nano-sized vesicles (j57| 
that have reported only prolate deformations. 

The oblate-prolate transition is independent of membrane properties; it is analogous to the one 
observed with drops (l58l |59|) . It is also independent of the viscosity ratio because the electrohydro- 
dynamic flow stops at steady state due the interface immobilization by gradients in the membrane 
tension. 

4.2.3 The effective dipole theory does not predict the prolate-oblate transition 



The effective dipole theory, summarized in Appendix C.2[ models the cell as a sphere with effective 
permittivity. The theory successfully explains the dielectrophoresis and electrorotation of cells, 
because it correctly describes the perturbation due to the cell in the exterior electric field. However, 
the internal electric field is not physical, which leads to incorrect interior Maxwell stress and 
electric force distribution on the membrane. Accordingly, the predicted deformation is oblate at 
low frequencies ()23p . which is at odds with the experimental observations with vesicles ([1]). 

Figure |9] compares the predictions of our model and the effective dipole theory for the electric 
tractions. It shows that the two models agree at frequencies u; > wi, where the transmembrane 
potential has vanished. At low frequencies, where the field inside the vesicle is zero, the effective 
dipole theory would correctly predict the electric tractions if only the contribution from the exterior 
electric field is taken into account. However, at intermediate frequencies, where the vesicle interior 
participates in the conduction process and the transmembrane potential is still significant, i.e., 
LO ^ uji^ the effective dipole theory diverges from our model as well as experimental observations ([T]). 
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5 Conclusions and outlook 

We have developed a theory that explains the observed morphological transitions of vesicles in 
a uniform AC electric field, in particular, the shape dependence on the field frequency and con- 
ductivity ratio between the inner and outer fluids. Prolate deformations at low frequencies have 
purely dielectric origin and result from electric pressure due to polarization charges pulling the 
vesicle at the poles. Oblate deformations, however, result from induced free surface charges, which 
cause negative pressure and transient electrohydrodynamic flow driven by tangential electric trac- 
tions. The prolate-oblate transition at low frequencies depends on the membrane capacitance and 
conductance. At high frequencies, electric stresses become negligible and do not aflFect the vesicle 
equilibrium quasi-spherical shape. The theory also predicts a high-frequency oblate-prolate tran- 
sition, which is analogous to the one observed with drops: it is independent of the membrane 
electric properties and depends only on the conductivities of inner and outer fluids. The transition 
frequency, however, is not given by the Maxwell- Wagner polarization time, but is determined by 
electrohydrodynamics . 

We have considered the problem for vesicle electrodeformation from a mechanical point of view 

where the vesicle shape is determined by the balance of forces exerted on the interface. Thus, our 
formalism can be easily extended to electric fields of arbitrary symmetry as well as to situations 
when external electric and fiow fields are simultaneously applied. 

Our current theory is a step in a systematic study of the electrohydrodynamics of deformable 
cells and, as such, some potentially important effects arc neglected. First, our treatment assumes 
that all media are electrically homogeneous and is based on solutions of Laplace's equation. This 
approach requires that the Debye length of the media is small compared to the radius of the vesicle 
or the thickness of the membrane. Thus our theory might break down at low conductivities and 
frequencies. Second, the model does not include shear elasticity of the membrane, which is essential 
in the mechanics of the red blood cell. Third, the membrane is assumed to be non-permeable to 
ions. However, at low frequencies, the duration of application of the electric field may be sufficient 
to porate the membrane. An electric current due to ion movement through field-induced pores 
would affect the electric field and tractions, and therefore vesicle shapes. Electrokinetic effects, 
the role of shear elasticity, membrane poration, and membrane charge represent interesting and 
challenging problems to be investigated in the future. 
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A List of symbols 

subscript r denotes radial 
subscript 9 denotes tangential 
sub/superscript "el" denotes electric 
sub/superscript "hd" denotes hydrodynamic 
sub/superscript "mm" denotes membrane 
sub/superscript "in" denotes interior 
sub/superscript "ex" denotes exterior 
superscript * denotes complex conjugate 
ReW denotes Real part of [ ] 
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/m[] denotes Imaginary part of [ ] 

a vesicle radius 
Cm membrane capacitance 
Ca capillary number 
E electric field 

fjm shape deformation parameter 

Gm membrane conductivity 

h membrane thickness 

H mean curvature 

p pressure 

S permittivity ratio 

t time 

T bulk stress 

V fluid velocity 

Vm transmembrane potential 

X = h/a dimensionless membrane thickness 

Yjm spherical harmonic 

T] viscosity 

p density A conductivity 
e permittivity 
A conductivity ratio 
X viscosity ratio 

Xs membrane viscosity parameter 

LOi frequency of the prolate-oblate transition 

UJ2 frequency corresponding to the inverse Maxwell- Wagner polarization time 

frequency of the oblate-prolate transition 
$ electric potential 
r tractions 
a membrane tension 
K bending rigidity 
A excess area 



B Spherical harmonics 

The normalized spherical scalar harmonics are defined as (60) 



Yjm {0, ip) 



2j+l (i-m)! 
47r (jr'+m)! 



(-l)'"P'"(cos6l)e''"'^, 



(47) 
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C Electrostatic field and stresses for a spherical shell 

C.l Our model: A sphere with interfacial capacitance and conductivity 

Schwann et al. (|4H [53|) have solved the problem for the electric field about a spherical shell with 
radius a and shell thickness h to obtain Eq. [8] for the potential difference between the inner and 
outer shell surfaces. Assuming a very thin shell h/a <^ 1, we can approximate the membrane 
with a two-dimensional interface that possesses capacitance. Accordingly, the spherical shell is 
approximated by a sphere with a discontinuous potential at the interface. 

Solving Eq. [7] with the boundary conditions Eq. [8j Eq. [9] and Eq. [TO] leads to 

= -[r + P<=V-2]exp(ia;t)cos6', = -P'"r exp(icjt) cos 6* . (49) 

where 

pex _ ^oo + fccx) + fejnKra pin _ jooj, 3 — 2Vm , . 

"'in ■^"'cx "'in ■^"'ex 

and k denote the dimensionless complex conductivities of the inner and outer fluids 

fcin = A + IUjS , fcex = I + liO . (51) 

The tractions are computed from the Maxwell stress tensor. The radial (pressure) component 
is given by 

<i = ^[-2{Tf + r|)5 + 5r| - 2d^T, + 5r| + 2{d^)^] , (52) 
and the tangential (shearing) component is 

= -i^ [{r? + r|)5 + 2r| + (d-)r3 + 2r| - (d-)^] , (53) 

where ri = i2e[P™] ,r2 = /m[P™],r3 = i?e[i-"^^] , T4 = Im[P'^^]. Re[] and /m[] denote real and 
imaginary part. Taking the zero-thickness limit, x = 0, our solution reduces to the result for a 
spherical drop ()58|) 

^ei,drop ^3(^i^^2_2S + {S- l)^Suj^) (2 + S)\uj^ + ujI)-^ , (54) 
^el4rop ^ _ ^)(2 + Snu;^ + u;i)-i , (55) 

where 102 is given by Eq. [37] 

The effective charge density is calculated from Eq. |9] 

Q{io, t) = Qci,'^) cos(wt) -|- (75(0;) sm{ut) , (56) 

where 

Qs(w) = 2t4 + St2 , qs{oj) =d^- 2t3 - Sn (57) 
The frequency dependence of the charge density can be cast into the form 

Q{LO,t) = Q{uj)cos{Lot + ip) , (58) 
where the amplitude is Q{uj) = [qs + Qc]"^ , and the phase shift is ip = Qs/qc- 
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C.2 The effective sphere model 

The dipole theory models the cell as a sphere with an effective permittivity (|13| [1 

-1 



eff 



(l-x)-3 + 2 



kin ~l~ '^kmm 



(l-x) 



-3 ^in ^mm 
kin ~l~ 2fejnm 



where 



Mn = A + iu;5 , A:mm = Amm + iwS'n 

The electric field is described by an electric potential 



r + f 



-2 ( ^il^ ~^ ^cx) 



ex J 

/Cj + 2A;ex 



COS I 



(59) 
(60) 

(61) 



The electric potential is continuous, and hence there is no transmembrane potential within the 
framework of the effective dipole theory. 



D Deformation of a prolate vesicle in strong fields 

Consider an initially non-spherical, non-fluctuating vesicle. This situation can occur in strong elec- 



tric fields, where the vesicle is already maximally deformed, Eq. 42 and then the field direction is 



changed. The evolution to the new stationary shape is no longer described by Eq. 44 because the 



tension is no longer given by Eq. 40 The effective tension has to be determined self-consistently 
along with the field-induced changes in shape to keep the total area constant (154]) , see Supple- 
mentary material. The leading order vesicle electrohydrodynamics becomes non-linear in contrast 
to the corresponding results for drops and capsules (l38l l58l |6T]) . This feature of non-equilibrium 
vesicle dynamics has been noted by several authors in relation to vesicle dynamics in shear flow 

The vesicle deformation described by Eq. l39] can be approximated by 



/20 = C^'(l-2A-Vlo) 



2m 



-2C7'='A-V2o/; 



2m 



(62) 



where the dot denotes time derivative. The modes /2m are slaved to the /20, which is forced to 



change by the electric field. Eq. 62 can be integrated to yield 



/2o(i) = (^tanh 



\ 



(63) 



where 5 is the maximum possible deformation 



17 



References 

1. Aranda, S., K. A. Riske, R. Lipowsky, and R. Dimova, 2008. Morphological transitions of 
vesicles induced by AC electric fields. Biophys. J. 95:L19-L21. 

2. Zhao, M., B. Song, J. Pu, T. Wada, B. Reid, G. R Tai, F. Wang, A. H. Guo, P. Walczysko, 
Y. Gu, T. Sasaki, A. Suzuki, J. V. Forrester, H. R. Bourne, P. N. Devreotes, C. D. McCaig, and 
J. M. Penninger, 2006. Electrical signals control wound healing through phosphatidylinositol- 
3-OH kinase- and PTEN. Nature 442:457-460. 

3. Funk, R. H. W., and T. K. Monsees, 2006. Effects of Electromagnetic Fields on Cells: Physio- 
logical and Therapeutical Approaches and Molecular Mechanisms of Interaction. Cells Tissues 
Organs 182:59-78. 

4. Voldman, J., 2006. Electrical Forces For Microscale Cell Manipulation. Annu. Rev. Biomed. 
Eng. 8:425-454. 

5. Zimmermann, U., and G. A. Neil, 1996. Electromanipulation of cells. CRC Press, Boca Raton. 

6. Neumann, E., A. E. Sowers, and C. A. Jordan, 1989. Electroporation and electrofusion in cell 
biology. Plenum Press, New York. 

7. Chizmadzhev, Y. A., P. Kuzmin, and V. P. Pastushenko, 1985. Theory of the dielectrophoresis 
of vescicles and cells. Biol. Mem. 2:1147-1161. 

8. Chizmadzhev, Y. A., P. Kuzmin, and V. P. Pastushenko, 1988. Dielectrophoresis and elec- 
trorotation of cells: unified theory for spherically symmetric cells with arbitrary structure of 
membrane. Biol. Mem. 5:65-78. 

9. Turcu, I., and C. M. Lucaciu, 1989. Dielectrophosresis -a spherical shell model. J. Phys. A 

22:985-993. 

10. Jones, T. B., 1995. Electromechanics of particles. Cambridge University Press, New York. 

11. Gimsa, J., and D. Wachner, 1999. A polarization model overcoming the geometric restrictions 
of the laplace solution for spheroidal cells: Obtaining new equations for field-induced forces 
and transmembrane potential. Biophys. J. 77:1316-1326. 

12. Dolinsky, Y., and T. Elperin, 2006. Dynamics of a spheroidal particle in a leaky dielectric 
medium in an ac electric field. Phys. Rev. E 73:066607. 

13. Foster, K. R., F. A. Sauer, and H. P. Schwan, 1992. Electrorotation and levitation of cells and 
colloidal particles. Biophys. J. 63:180-190. 

14. Jones, T. B., 2003. Basic theory of dielectrophoresis and electrorotation. IEEE Eng Med Biol 
Mag. 22:33-42. 

15. Miller, R. D., and T. B. Jones, 1993. Electro-orientation of ellipsoidal erythrocytes. Theory 

and experiment. Biophys. J. 64:1588-1595. 

16. Engelhardt, H., and E. Sackmann, 1988. On the measurement of shear elastic moduli and 
viscosities of erythrocyte plasma membranes by transient deformation in high frequency electric 
fields. Biophys J. 54:495-508. 



18 



17. Bryant, G., and J. Wolfe, 1987. Electromechanical stresses produced in the plasma membranes 
of suspended cells by applied electric fields. J. Membr. Biol. 96:129-139. 

18. Poznanski, J., P. Pawlowski, and M. Fikus, 1992. Bioelectrorheological model of the cell. 3. 
Viscoelastic shear deformation of the membrane. Biophys. J. 61:612-620. 

19. Helfrich, W., 1974. Deformation of lipid bilayer spheres by electric fields. Z. Naturforsch. 
29c:182-183. 

20. Kunimrow, M., and W. Helfrich, 1991. Deformation of giant lipid vesicles by electric fields. 
Phys. Rev. A 44:8356-8360. 

21. Winterhalter, M., and W. Helfrich, 1988. Deformation of spherical vesicles by electric fields. J. 
Coll. Int. Sci. 122:583-586. 

22. Teissie, J., and T. Y. Tsong, 1981. Electric field induced transient pores in phospholipid bilayer 
vesicles. Biochemistry 20:1548-1554. 

23. Sukhorukov, V., H. Mussauer, and U. Zimmermann, 1998. The effect of electrical deformation 
forces on the electropermeabilization of erythrocyte membranes in low- and high-conductivity 
media. J. Memb. Biol. 163:235-245. 

24. Isambert, H., 1998. Understanding the electroporation of cells and artificial bilayer membranes. 
Phys. Rev. Lett. 80:3404-3407. 

25. Sens, P., and H. Isambert, 2002. Undulation instability of lipid membranes under an electric 
field. Phys. Rev. Lett. 88:Art. No. 128102. 

26. Lacoste, D., M. Lagomarsino, and J. Joanny, 2007. Fluctuations of a driven membrane in an 
electrolyte. Europhys. Lett. 77:18006. 

27. Ambjornsson, T., M. A. Lomholt, and P. L. Hansen, 2007. Applying a potential across a 
biomembrane: Electrostatic contribution to the bending rigidity and membrane instability. 
Phys. Rev. E 75:051916. 

28. Krassowska, W., and P. Filev, 2007. Modeling electroporation in a single cell. Biophys. J. 
92:404-417. 

29. Riske, K. A., and R. Dimova, 2005. Electro-deformation and poration of giant vesicles viewed 
with high temporal resolution. Biophys. J. 88:1143-1155. 

30. Riske, K. A., and R. Dimova, 2006. Electric pulses induce cylindrical deformations on giant 
vesicles in salt solutions. Biophys. J. 91:1778-1786. 

31. Dimova, R., S. Aranda, N. Bezlyepkina, V. Nikolov, K. A. Riske, and R. Lipowsky, 2006. A 
practical guide to giant vesicles. Probing the membrane nanoregime via optical microscopy. J. 
Phys. Cond. Matt. 18:S1151-S1176. 

32. Dimova, R., K. A. Riske, S. Aranda, N. Bezlyepkina, R. L. Knorr, and R. Lipowsky, 2007. 
Giant vesicles in electric fields. Soft matter 3:817-827. 

33. Mitov, M. D., P. Meleard, M. Winterhalter, M. I. Angelova, and P. Bothorel, 1993. Electric- 
field-dependent thermal fluctuations of giant vesicles. Phys. Rev. E 48:628-631. 



19 



34. Peterlin, P., S. Svetina, and B. Zeks, 2007. The prolatc-to-oblatc shape transition of phospho- 
hpid vesicles in response to frequency variation of an AC electric field can be explained by the 
dielectric anisotropy of a phospholipid bilayer. J. Phys. Cond. Phys. 19:136220. 

35. Hyuga, H., K. Kinosita Jr., and N. Wakabayashi, 1991. Deformation of vesicles under the 
influence of strong electric fields. Jpn. J. Appl.Phys. 30:1141-1148. 

36. Nelson, P., 2004. Life in the Slow Lane: The Low Reynolds-Number World. In Biological 

Physics: Energy, Information, Life, Freeman, 158-194. 

37. Sauer, F. A., 1985. Interaction forces between microscopic particles in an external electro- 
magnetic field. In A. Chiabrera, C. Nicolini, and H. P. Schwan, editors, Interactions between 
electromagnetic fields and cells. Plenum Press, 181-202. 

38. Taylor, G. I., 1966. Studies in electrohydrodynamics. I. Circulation produced in a drop by an 
electric field. Proc. Royal Soc. A 291:159-166. 

39. Melcher, J. R., and G. I. Taylor, 1969. Electrohydrodynamics - a review of role of interfacial 
shear stress. Annu. Rev. Fluid Mech. 1:111-146. 

40. Saville, D. A., 1997. Electrohydrodynamics: The Taylor-Melcher leaky dielectric model. Annu. 
Rev. Fluid Mech. 29:27-64. 

41. Schwan, H. P., 1989. Dielectrophoresis and rotation of cells. In E. Neumann, A. E. Sowers, and 
C. A. Jordan, editors, Electroporation and electrofusion in cell biology. Plenum Press, 3-21. 

42. Kinosita Jr., K., I. Ashikawa, N. Saita, H. Yoshimura, H. Itoh, K. Nagayama, and A. Ikegami, 
1988. Electroporation of cell membrane visualized under a pulsed laser fluorescence microscope. 
Biophys. J. 53:1015-1019. 

43. Needham, D., and R. M. Hochmuth, 1989. Electromechanical permeabilization of lipid vesicles. 
Role of membrane tension and compressibility. Biophys. J. 55:1001-1009. 

44. Leal, L. G., 1980. Particle motions in a viscous fluid. Ann. Rev. Fluid Mech 12:435-476. 

45. Sozou, C, 1972. Electrohydrodynamics of a liquid drop - time-dependent problem. Proc. Royal 
Soc. A 331:263-272. 

46. Leal, L. G., 1992. Laminar Flow and Convective Transport Processes. Butterworth-Heinemann, 
Boston. 

47. Petrov, A. G., 2006. Electricity and mechanics of biomembrane systems: Flexoelectricity in 
living membranes. Analytica Chimica Acta 568:70-83. 

48. Raphael, R. M., A. S. Popel, and W. E. BrowneU, 2000. A Membrane Bending Model of Outer 
Hair Cell Electromotility. Biophys. J. 78:2844-2862. 

49. Seifert, U., 1999. Fluid membranes in hydrodynamic flow fields: Formalism and an application 
to fluctuating quasispherical vesicles. Eur. Phys. J. B 8:405-415. 

50. Seifert, U., 1997. Configurations of fiuid membranes and vesicles. Advances in physics 46:13- 
137. 

51. Discher, D. E., and F. Ahmed, 2006. Polymersomes. Annu. Rev. Biomed. Eng. 8:323-341. 



20 



52. Dimova, R., U. Seifert, B. Poligny, S. Forstcr, and H.-G. Dobereiner, 2002. Hyperviscous 
diblock coploymer vesicles. Eur. Phys. J. D 7:241-250. 

53. Grossc, C, and H. P. Schwan, 1992. Cellular membrane potentials induced by alternating 
fields. Biophys. J. 63:1632-1642. 

54. Vlahovska, P. M., and R. Gracia, 2007. Dynamics of a viscous vesicle in linear flows. Phys. 
Rev. E 75:016313. 

55. Evans, E., and W. Rawicz, 1990. Entropy driven tension and bending elasticity in condensed- 
fluid membranes. Phys. Rev. Lett. 64:2094-2097. 

56. Dimova, R., B. Poligny, and C. Dietrich, 2000. Prctransitional effects in DMPC-vesicle mem- 
branes: optical dynamometry study. Biophys. J. 79:340-356. 

57. Kakorin, S., T. Licsc, and E. Neumann, 2003. Membrane curvature and high-fleld electropora- 
tion of lipid bilayer vesicles. J. Phys. Chem.B 107:10243-10251. 

58. Torza, S., R. Cox, and S. Mason, 1971. Electroliydrodynamic deformation and burst of liquid 
drops. Phil. Trans. Royal Soc. A 269:295-319. 

59. Vizika, O., and D. A. Saville, 1992. The electrohydrodynamic deformation of drops suspended 
in liquids in steady and oscillatory electric fields. J. Fluid Mech. 239:1-21. 

60. Varshalovich, D. A., A. N. Moskalev, and V. K. Kheronskii, 1988. Quantum Theory of Angular 
Momentum. World Scientfic, Singapore. 

61. Ha, J. W., and S. M. Yang, 2000. Electrohydrodynamic effects on the deformation and orien- 
tation of a liquid capsule in a linear flow. Phys. Fluids 12:1671-1684. 

62. Misbah, C, 2006. Vacillating breathing and tumbling of vesicles under shear flow. Phys. Rev. 
Lett. 96:028104. 

63. Lebedev, V. V., K. S. Turitsyn, and S. S. Vergeles, 2007. Dynamics of nearly spherical vesicles 
in an external flow. Phys. Rev. Lett. 99:218101. 



21 



FIGURE LEGENDS 



Figure [TJ An illustration of the streamlines of the electrohydrodynamic flow inside a drop, sur- 
face charge distribution. The corresponding direction of the tangential electric traction is denoted 
by arrows (I38| l58]l . (a) interior fluid less conducting than the exterior one, A/S < 1; (b) interior 
fluid more conducting than the exterior one, A/S > 1. 

Figure [2j A sketch of a vesicle in a uniform electric field. The zoomed region of the interface 
illustrates the lipid bilayer structure of the membrane. 

Figure [sj (a) The transmembrane potential at the poles = 0, vr calculated from Eq. 32 



and (b) the interior electric field calculated from Eq. 50 for a spherical shell with membrane 
thickness x = 5 x 10~^ (solid line) and x = 10~^ (dashed line) in a uniform AC electric field. 
A = 0.5, S = 1.001, Cm = 0.025/x, Gm = 0. 

Figure |4} Electric tractions as a function of frequency for conductivity ratios A = 0.5 (solid 
line) and A = 1.5 (dashed line). The other parameters are S = 1.001, Cm = 0.025/x, Gm = and 
membrane thickness x = 5 x 10~^. Dotted lines represent the electric tractions in the case of a 
droplet (the zero-membrane-thickness limit). 



(a) Electric pressure calculated from Eq. 52 



(b) Tangential electric force calculated from Eq. 53 



Figure [5j The absolute value of the amplitude of the interfacial charge density at the poles 
9 = 0,11. Parameters are as in Fig. |4j only A = 0.5 is plotted. The dashed line represents the 
induced charge in the case of a droplet (the zero-membrane-thickness limit). 

Figure |6j Evolution of the ellipsoidal deformation /20, calculated from Eq. 44 of a quasi- 



spherical vesicle upon application of a uniform AC electric field. Parameter values are A = 0.2, 
A = 1.5, X = 5 X 10~'^, Gm = 0, and Gm = 50. Ca = 684, which corresponds to Eq = W^V/m, 
rjcx = 10~^Pa.s, K = 25kBT and a = lOfim. 

(a) for frequencies uj = 0.01, 0.1, 1, denoted by solid, dashed and long-dashed lines, viscosity ratio 
X = 1 and Xs = 0. 

(b) for viscosity ratios x = Oi 5, denoted by solid, dashed and long-dashed lines, at AC field 
frequency uj = 0.01 and Xs = 0. 

(c) for membrane viscosity parameters Xs = 0, 10, 100 denoted by solid, dashed, and long-dashed 
lines, at AC field frequency u = 0.01 and viscosity ratio x = 1- 

Figure [7j Ellipsoidal deformation /2o//2o'^^ vesicles in AC field as a function of frequency. 
Solid line is for conductivity ratio A = 0.5 and dashed line is for A = 1.5. Other parameters are 
A = 0.2, X = 1, Xs = 0, X = 5 X 10-^ Gm = and Gm = 50. 

Figure [8j Prolate-oblate transition frequency for different conductivity ratios A. The points 
are experimental data ([ll [32|l for interior fluid conductivities of about lO^^S/m averaged over 
number of vesicles with different size. The solid line is the theoretical prediction using x = 10~^ 
corresponding to vesicle radius a = 50/im. The bottom dashed line is calculated using x = 2 x 10~^ 
and the top dashed line is calculated using x = 5 x 10~^ . The basic charging time is estimated to 
be tc = W~'^s. 

Figure [9] The electric tractions according to our model (solid line), the drop model (short- 
dashed line), the effective dipole theory (long-dashed line) and only exterior field contribution in the 
effective dipole theory (dot-dashed line), (a) electric pressure (b) tangential tractions. Parameters 
are: A = 0.5, A = 1.5, Xs = 0, x = 5 x 10-^ Gm = 0,Cm = 50 and Ca = 684. 
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